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Studying a certain sub class of higher order Horndeski (scalar-tensor) theories we discuss a method 
discovered recently permitting analytic black hole solutions with a non trivial and regular scalar field. 
One of the solutions found has de Sitter asymptotics and self tunes the bulk cosmological constant. 
Using the aforementioned method we find and analyse new black hole solutions which can also have 
de Sitter asymptotics. By looking at small deviations of the integration constant responsible for self 
tuning we discuss the robustness of the self tuning mechanism. We find that neighboring solutions 
to the one previously found present also self tuning properties with unaltered effective cosmological 
constant. 


I. INTRODUCTION-MOTIVATION 

General Relativity (GR) is a theoretically and experimentally consistent theory. It is also mathematically unipue. 
Indeed theoretical consistency and mathematical uniqueness emanate as a corollary to Lovelock’s theorem^ [l| in 
4 dimensions: consider a smooth Lagrangian depending up to second order derivatives of the metric tensor C = 
C{M., g ,V g ,W g) and endowed with a Levi Givita connection. Then in D = 4 dimensions GR with cosmological 
constant is the unique metric theory emerging from the action, 

%)=/ d^x^/^)[R-2K\ ( 1 . 1 ) 

Jm 


which upon metric variation yields, 

• Equations of motion which are 2'^'^-order PDE’s, 

• given by a symmetric two-tensor, G^i, + Kg^i, 

• and admitting Bianchi identities. 

Under these hypotheses GR with a cosmological constant is the unique massless, spin 2, 4 dimensional theory of 
gravity! In other words the Einstein-Hilbert action is the unique metric action which evades Ostrogradski no-go 
theorem of 1850 j^. This no-go theorem states that any non-degenerate, higher than second order derivative theories 
inevitably lead to ghost degrees of freedom! In fact the only safe metric term we can add to the above action 
is the Gauss-Bonnet scalar, E? — ARiccP + Riemann^, whose action is a 4 dimensional topological invariant the 
generalized Euler characteristic. Although this geometric term is not trivial its metric variation is not modifying thus 
GR equations of motion. 

GR is a physical theory which has remarkable agreement with weak and strong gravity experiments from scales 
smaller than the millimeter and up to 30 times Earth to Sun distance... At astrophysical scales however, GR equations 
of motion, require dark matter interacting only gravitationally. The presence of dark matter although well founded 
theoretically still evades astroparticle and accelerator experimental detection. Furthermore, at cosmological scales 
the Universe is now observed to be in an accelerating phase requiring some source of yet unknown dark energy in 
order for Einstein’s equations to agree with a multitude of differing cosmological and astrophysical observations 
In resume, if we assume only ordinary (detected) sources of matter there is a flagrant disagreement between local, 
astrophysical and cosmological data. In fact roughly about 4% is matter of known source in the Universe. Even if we 
accept theoretically the presence of dark matter, dark energy of about 70% still has to be accounted for. For dark 


*Electronic address: christos.charmousis@th.u-psud.fr 
^Electronic address: diosifid@auth.gr 

^ The four dimensional version of this theorem was already known early on Q] and mathematically stems from the mathematical works 
of Chern 0 




2 


energy, which we will focus from here onwards, there is an almost trivial way out. Assume a cosmological constant of 
energy density p\ = = (10“^ey)^, ie. modify Einstein’s equation by, 

G ^1/ ^obsSfiu — ( 1 -^) 

Now this is a rather economic and consistent solution since a cosmological constant is the only other metric term 
allowed given the above uniqueness theorem dm. However this cosmological scale is tiny, the mass scale corresponds 
to neutrino mass scales. This is because a cosmological constant introduces a scale in the metric, tq ^ ^/A which is 
as tiny, as the inverse size of the Universe today, tq = Note that the quotient made out of Solar system scales 

divided by Cosmological Scales leads to an enormous number = 10“^® at which we are probing the physical 

theory in question, GR, compared to the local experimentally robust scales. In fact suppose we were to use local 
experiments in order to bound the cosmological constant value (without taking into account cosmological data). That 
would involve constraining the cosmological constant from the (very) well known trajectories of planets. In essence, 
replacing Schwarzschild geometry by de Sitter Schwarzschild and taking the PPN limit. Local experiments would 
allow a cosmological constant roughly I0^° bigger than its current cosmological bound! Physical theories describe the 
physical world from good to remarkable precision but up to (or at) certain scales. Are we therefore beyond the scale 
of application of GR? Although deep in the UV, i.e. close to the Planck scale, GR is ill defined, could it be that also 
at vast and extremely dilute cosmological scales, of almost zero curvature, some other theory completes the theory of 
General Relativity? Furthermore, the cosmological constant is accompanied by the cosmological constant problem: 
its huge theoretical value, as expected from QFT considerations is at complete contrast with a finely tuned and tiny 
cosmological constant @ (for a complete discussion and interesting developments see also i). 

Summarizing: what if the need for exotic matter or a tiny cosmological constant is the sign for novel gravitational 
physics at very low energy scales or large distances. It is clearly a valid avenue of theoretical research in order to 
understand the robustness of GR at 10^® greater than local distances. In fact historically there is a very similar 
paradigm encountered at the advent of GR. More than 100 years ago, there were also observational short comings of 
Newton’s theory of gravity. Indeed, Mercury the closest orbiting planet around the sun was deviating from Keplerian 
motion by a tiny yet definite advance in its perihelion. At the time, differing possibilities were pursued including most 
interestingly, the existence of a novel planet, Vulcan, which was postulated to orbit even closer to the sun, slightly 
deviating the perihelion of Mercury much like a dark source to Newtonnian equations. In fact a simple next order 
correction with one additional parameter (much like a cosmological constant) was enough to save Newton’s laws (at 
the experimental precision of the time). Furthermore, although the presence of Vulcan may seem now days an ad 
hoc hypothesis it actually was not. The outermost planet of the solar system, Neptune, was in fact first theoretically 
predicted prior to being observed by the (same) astronomer Le Verrier. Indeed the theoretical prediction was triggered 
due to the reported inconsistencies in the Newtonian planetary trajectory of Uranus-which were explained by the 
presence of Neptune (see the exciting account in 0)1 For the advance of the perihelion of Mercury, however, it 
was not the ’’simpler” solution that was the correct one, indeed GR, a far more involved than Newtonian theory 
had to be employed in order to explain strong gravity at local scales as was accepted by the scientihc community 
only later by Eddington s famous eclipse experiment. Just as the eventual success of GR was not explaining the 
advance of Mercury’s perihelion, modification of gravity cannot only be ”an explanation” of the effective cosmological 
constant. Successes of GR underlies in all the theoretical and subsequently observational advances such as, the Big 
Bang Universe, inflation theory, black holes, binary pulsars, quasars etc. 

So which avenues should one choose in an effort to modify such a well established theory, both theoretically and 
experimentally? For a start, we can deduce a set of basic facts/guidelines that can serve as a rough guide to possible 
research directions: 

• Since GR is the unique 4 dimensional metric theory we will need to introduce new and genuine gravitational 
degrees of freedom! 

• Novel and non degenerate gravitational degrees of freedom must not lead to higher derivative equations of 
motion. For then additional degrees of freedom are ghosts and vacuum is unstable due to Ostrogradski theorem 
1850 as explained in (Tlj |. 

• Matter must not directly couple to novel gravity degrees of freedom. Matter sees only the metric and evolves 
in metric geodesics. As such EEP is preserved and space-time can be put locally in an inertial frame. 

• Novel degrees of freedom need to be screened from local gravity experiments. We need a well defined GR local 
limit such as the Ghameleon, see for examplefl^, or the Vainshtein mechanism M- The equivalent of this in 
GR would be the weak gravity Newtonian limit. 
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• Exact gravitational solutions are essential in modified gravity theories in order to understand strong gravity 
regimes and novel characteristics. We will need to deal with the no hair paradigm and the absence of GR 
theorems such as Birkhoff’s theorem in modified gravity theories. 

• A modified gravity theory should tell us something about the short-comings of GR such as the cosmological 
constant problem and in particular how to screen an a priori enormous cosmological constant. 

The above list of points are not exhaustive and one must treat them not dogmatically but rather as rough guidelines. 
For example f(R) theory has higher order derivative equations of motion but does not necessarily acquire ghost degrees 
of freedom. The theory can be easily transformed in a scalar-tensor theory. The spin 2 degrees of freedom are in fact 
not present for this (degenerate) case but only an extra scalar degree of freedom. Having this in mind we can sketch, 
again with precaution, differing directions to go beyond GR (for a full account see for example 0 ): 

• We can assume metric theories but in the presence of extra dimensions : Extension of GR to Lovelock theory 
with modified yet second order field equations [l3|. Here belong braneworld models such as the DGP model 
M , RS models [l^ or the generic case of Kaluza-Klein compactification. 

• One can assume that the graviton is not massless but massive. In fact in this case generically graviton potentials 
suffer from a Boulware Deser ghost mode [i3- The only possible theories avoiding this specific mode are non 
trivial and are parametrized by 3 independent couplings (see for example 0 ). Hence massive gravity and 
bigravity theories are greatly constrained already theoretically and the most general interaction (potential) 
term is known. 

• One can assume a 4-dimensional scalar modification of GR. This is the case of Scalar-tensor theories fl^ . [2lj, 
f{R) [1^, Galileon/Hornedski theories [l^, [2^,[2^ (for a review and relation to Lovelock theory see |26j|'). 

• We can consider breaking symmetries present in GR such as Lorentz invariance symmetry. This has interesting 
implications in the UV as it is a way of obtaining higher order theories which are power counting renormalizable 
and have a priori better ultra-violet properties: this is the case of Horava gravity [l^, [1^ and in its infra red 
limit, Einstein Aether theories [2^, which are healthy and unique vector tensor theories with a consistent GR 
limit. 

• Theories modifying geometric properties of spacetime such as with the inclusion of torsion or a different choice 
of geometric connection (see for example [13 )■ 

In this paper we will consider the simpler of modifications with a single extra degree of freedom, scalar-tensor theory. 
Scalar tensor theories are in a certain sense a prototype of gravity modification and have some nice properties. In 
particular, 

• They admit a uniqueness theorem due to Horndeski, 1973, which much like Lovelock’s theorem determines the 
most general form of their action [13, [23| . 

• They contain or are limits of other modified gravity theories. f{R) theory for example is a scalar tensor theory 
in disguise whereas massive gravity or bigravity are scalar tensor theories in their decoupling limit. This is also 
true of certain higher dimensional braneworlds such as DGP and even for Kaluza-Klein compactifications of 
Loveock theory [^ , [13 , [13 

• They can have insightful screening mechanisms as in the Chameleon [13 and Vainshtein [l3 mechanisms. 

• They include specific Galileon terms, dubbed Fab 4, that can screen classically a big cosmological constant [13. 

• They specifically break GR black hole uniqueness theorems such as Birkhoff’s theorem and suffer from no-hair 
black hole theorems. 

It is this latter characteristic we will study and try to overcome for we will seek black hole solutions of such theories 
where the scalar field will not be trivial but well behaved. Before doing so let us point out the most general theory 
at hand. Horndeski showed several years ago with techniques similar to those of Lovelock that the most general 
scalar-tensor action with second order field equations is given by 


Sh 


-I 


(R^Xy/—g (L 2 -|- L 3 -|- L 4 + L 5 ) 
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L2 = K{cl>,X), 

Ls = -G^{ct>,X)V^, 

U = Gi{<P, X)R + G4x [(v2<^)2 - [V^VAf] , 

Lg = G5(0, ^ [(V2,^)3 - 3V2,^(V^V,</.)" + 2(V;,V,.^)3] 

where the Gi are unspecified functions of (p and X = and Gix = dGijdX. In fact this is the action of 

covariant Galileons found independently in [s^ in a different and elegant manner. Due to the higher order kinetic 
terms in X the theory generically screens the scalar mode locally following the Vainshtein mechanism (T^ . 

In the next section we will discuss scalar tensor black holes starting with a brief discussion of the no hair paradigm, 
and then certain solutions some of which are new and appear in the last subsection. We will then briefly conclude. 


II. SCALAR TENSOR BLACK HOLES 

Modified gravity theories have to deal with the absence of black hole uniqueness theorems (present in GR) as well as 
no hair theorems. Indeed, ’’Black holes have no hair”, is the celebrated phrase conjectured long ago by Wheeler. The 
no hair paradigm takes concrete form as a set of mathematical theorems governing asymptotically flat (or de Sitter) 
and stationary black holes under certain hypotheses. The conjecture itself states, that apart from charges measured 
at infinity by a far away observer, no additional degrees of freedom can describe the black hole geometry (for a review 
see [1^). The physical picture is roughly the following: as black holes are formed they expel or eat up surrounding 
matter ending up, in their stationary phase, as rather ’’simple” bald objects characterized by a finite number of 
measurable charges at infinity: mass, angular momentum, electric and magnetic charge. A black hole, is a rather 
blunt gravitational object having specific charges and not allowing additional parameters primary hair which are 
not associated to a conserved charge, or, secondary hair of no additional parameters but non-trivial fields interacting 
with the black hole spacetime. An ideal observer measuring these charges out at infinity ’’knows” specifically which 
black hole she or he is dealing with. They do not have to venture close up in order to gather details about the black 
hole structure, its horizon properties etc. This is true for asymptotically fiat and stationary black holes under certain 
hypotheses^. As such, any scalar field, part of a scalar tensor theory is expected to be frozen to a constant in order 
to evade a naked singularity. In other words scalar-tensor theories in their prototype Brans Dicke form [13 basically 
have the same GR black holes with the scalar field fixed to a constant value. As such a spherically symmetric BD 
black hole is again a Schwarzschild black hole with a constant scalar field. Otherwise the scalar tensor spherically 
symmetric solutions are singular geometries describing rather stars which in BD theory are different solutions to those 
describing black holes. We emphasize here the difference with GR where the Schwarzschild solution is the exterior 
metric for a star and a black hole due to Birkhoff’s uniqueness theorem. This is not true in scalar tensor theories and 
it is a source of problems of scalar tensor theories with local experiments. There is a non-trivial example of a black 
hole geometry that almost escapes this rule and that is the BBMB black hole which we will briefly turn to now. 


A. An example: The BBMB black hole 

Consider the following scalar tensor action, 

( 16 ^ ~ ~ 

The particular feature of this action is that the scalar (p couples conformally to the metric. Indeed there is an invariance 
of the p equation of motion (EoM) under the conformal transformation 

9aj3 ' ^ l/a/5 — ^ l/a/3 (^■^) 

p ^ p = Vl~^P (2-3) 

The whole action is not conformally invariant. In fact the Einstein-Hilbert term explicitly breaks conformal invariance. 
However, for this action for which there is a remnant of conformal symmetry there exists a vacuum black hole geometry 
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Recently an interesting counter-example to black hole no hair theorems was found for stationary black holes with a massive complex 
scalar field. There crucially the scalar field was assumed not to have the same symmetries as the spacetime metric |20|| . 
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with non-trivial scalar field and secondary black hole hair. This is the BBMB solution a static and spherically 
symmetric solution, 


A TON 
\ r ^ 


ds^ = - 1- df + 


dr^ 


{i-fY 


(d6»2 


Sin 


' dd(/7^) 


with a non trivial scalar field 


0 = 


dTrG r — : 


This particular geometry is also encountered in GR. It is that of an extremal Reisner-Nordtstrom solution with charge 
equal to mass. Our ideal observer sitting far out of this scalar tensor black hole would observe a static black hole of 
mass TO other than Schwarzschild but with no electric charge. Indeed note that the scalar decays at infinity and does 
not carry a specific scalar charge. Hence we are not dealing with primary hair (by dehnition) rather secondary hair 
as the scalar couples non-trivially (conformally) to the metric modifying its geometry. Indeed the only parameter of 
the solution is the black hole mass to. Almost immediately however we encounter the shortcoming of the solution: 
the scalar field explodes at the horizon location at r = to. This can be remedied in the presence of a cosmological 
consta nt [3^ and a (j)^ conformal potential. The black hole can be even in an accelerated phase or have a Taub NUT 
charge [^. Here, for the asymptotically flat case, the singular scalar at the horizon is a problem and it is not clear 
whether the solution is a black hole i.e. an endpoint of matter collapse. Going to the Einstein frame we can actually 
see that the solution hides additional singularities there at r = 2m where the conformal transformation in between 
the frames is singular. In fact at r = 2m the (j) dependent factor multiplying R in the action is exactly zero and 
for r < 2m changes sign hinting on gravity acquiring a ghost kinetic term! Could self gravitating matter go beyond 
the r = 2m surface? Or again, is the BBMB solution an endpoint of gravitational collapse? It would seem unlikely 
and this is an open problem which could be resolved numerically determining the exact nature of this solution. The 
BBMB metric is however a unique and simple 4 dimensional example with a non trivial scalar field. We will proceed 
now to remedy at least some of these problems with higher order Galileon fields. 


III. CONSTRUCTING GALILEON BLACK HOLES 


We will now consider the action^, 

S = [ d'^xy/^ (R- 2A — T] {dcj)Y + d^(j)du4> ; 


(3.1) 


which is part of the Horndeski action given above (HID comprising of the Einstein-Hibert term with a cosmological 
constant as well as the canonical scalar kinetic term and a higher order scalar tensor interaction term. The latter 
term is a higher order Galileon, central in our forthcoming discussion, yielding second order field equations. This nice 
property boils down to the divergence freedom of the Einstein tensor. Finally C, t] and /3 are coupling constants in 
order to keep track of the terms in the action. Indeed the metric field equations are. 


CG, 


fiiy 


-■q 


where is the double dual curvature tensor given by, 


Pfipua — Pppva PypSp’fT RpcrSi/p RpaQpv Rp,v9p<7 H” per Sp’^rSi^p^’ 


(3.2) 

(3.3) 

(3.4) 

(3.5) 


This tensor is divergence free has the same symmetry properties as the Riemann tensor and tracing two of its indices 
gives, 




= -G 


pa 


^ The details summarized in this section can be found in the original paper Q. See also 0 for the first black hole solutions discussed 
for this theory 
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the Einstein tensor. 

An important property of the above action is that the scalar field has translational invariance: (/) —>■ 0+const. This 
is important for then the scalar field equation can be written in terms of a conserved vector current, 

= 0, d^cj). (3.6) 

Note in the last expression of the current above, the explicit appearance of the metric field equations of the lower 
order Einstein-Hilbert and cosmological constant terms. Schematically we have that, R d^(j)dv(t> whereas 

A —^ d^4>dv<t’ associating thus in this intriguing way the terms present in the action. Summarizing therefore, 

apart from the latter property, the above action crucially includes higher order Galileon terms and has translational 
invariance for the scalar field. 

We will consider a static and spherically symmetric spacetime, 

ds^ = -h{r)dt^ + + r'^dV? (3.7) 

f[r) 

parametrized by the two radial functions / and g. It is then evident that cj) = (j){r) then the scalar equation is 
integrable and reads, 

(773- - PGnV94>' = c (3.8) 

where c the integration constant stands for primary hair for a scalar-tensor black hole. A non zero value of c comes 
along with regularity problems since the current norm is singular = J^J'^g^u = {J’^Y'grr unless J’' = 0 at the 
horizon and hence c = 0. It is then easy to note that generically (j) = constant everywhere as noted in and we 
have the appearance of a no-hair argument i.e. the presence of a trivial scalar and a GR black hole solution. However, 
here there is a slight twist in this argument. Indeed, due to the higher order nature of the terms in the action we can 
instead of (/> = constant require that, 

/3G- - 773 - = 0. (3.9) 

and thus search for a different branch of solutions. Here note that the presence of the Einstein tensor (13.91) in the 
higher order term, is just the Hamiltonian constraint for a GR plus cosmological constant theory, and this is the 
reason for our Gallileon choice. A second point is that we can also include a linear time dependence in the scalar field 
(j){t,r) = qt + without affecting the scalar field equation, —dr[{f3G^'" — r]g^'")dr'il}] — dt[{f3G** — r]g**)dt{qt)\ = 0 
which is still verified given (13.9p . In fact, we can show that choosing (f>{t,r) = qt + '4>{r) solves also the (tr)-equation 
of motion and this at the end tells us that our Anzatz is mathematically consistent. This interesting property where 
the scalar field does not have the same symmetries as the background metric while having no flux is very similar to 
the idea used in [ 2 ^ to construct stationary minimally coupled scalar hair. We must stress however that although 
the scalar is not necessarily constant we have no scalar charge i.e. c = 0. Furthermore we have a finite current. We 
have satisfied the no hair argument not in permitting a non trivial scalar charge in the face of the constant c, but 
still having a scalar field which is not trivial. This is due to the higher order nature of the Galileon term providing us 
with an alternative geometric resolution /3G''’' — 773 '"’^ = 0 to the scalar and tr-equation. Therefore the fr-equation as 
well as the scalar field equation boil down to the geometric constraint, 

(/3 -I- ■qr‘^)h 
(3{rh)' 

fixing the spherically symmetric gauge and ^(t, r) = qt -\- '(/)(r). 

This solves therefore for the time dependence of </> and for the gauge function /. 
now need to solve for 7 />(r) and h(r) with the last two independent ODE’s, the (rr) 

(rr)-component we get ip' 

with 


In order to find solutions we 
and {tt) equations. From the 


X= (t] + piA 


. We see that for 77 = A = 0 time dependence of the scalar is essential in order to have a non trivial asymptotically 
flat solution. Finally the (tt)-component gives h{r) via. 


h{r) 



k{r) 


dr, 


r 


r 


P + rjr^ 


(3.11) 
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with 


g2^(/3 + - {KP + (2C^ - A) r^) k + Cofc^/^ = 0, (3.12) 

Therefore any solution to the algebraic eq for k = k{r) gives us a full solution to the system. For 77 = A = 0 we can 
easily find asymptotically flat solutions with a scalar field since then k is simply a constant which can be gauged away 
in a redefinition of time. The geometry is that of a Schwarzschild black hole Q and the scalar field due to the linear 
time dependence is regular at the black hole future horizon. The scalar field linearly diverges at infinity but note that 
its derivatives appearing in the action are well behaved. The solution is therefore a valid and regular solution unlike 
the BBMB case where the scalar blows up at the horizon location. 

One can then seek asymptotically de-Sitter solutions @ setting f = h whereupon 


k{r) 


(f3 + 

P 


(3.13) 


with 


= {Cv + P^)/{Pv), Cq = {(rj -/3A)^/^/r]. (3.14) 

The line element is that of de Sitter Schwarzschild f = h = 1— ^ with 'll;' = ±^\/l — h and cj){t, r) = qt + 'ip(r). 

Going to Eddington-Filkestein coordinates we can check that the solution is regular at the horizon. Note again that 
the scalar itself diverges at the de Sitter horizon but not its derivatives. 

The solution above has an interesting property : it is a GR like solution but crucially it does not directly depend 
on the vacuum cosmological constant A but rather on an effective cosmological constant Ae// = —r]/(3. This is true 
as long as the integration constant q^rj = A — Ae/y > 0 which is just p.l4l) . Also note that the solution belongs to 
a different branch than the GR solution of de Sitter Schwarzschild, one cannot switch off the higher order Gallileon 
term. Hence any arbitrary bulk cosmological constant of the action A > Ag//, which can be as big as the Planck scale, 
fixes q our integration constant leaving us with an effective geometric acceleration Ag// which can be tuned once and 
for all in the theory to be small. In other words the scalar tensor solution self tunes the vacuum cosmological constant 
as long as A > Ag//. But what is the behavior of other solutions to the algebraic equation (13.1211 ? In particular is 
there a fine tuning in the integration constant q? If we deviate from the conditions (13.141) will we still have de Sitter 
Schwarzschild asymptotics? We will turn to this question now. 


IV. GENERAL SOLUTIONS 

In the previous section we summarized the method and certain key solutions of spherical symmetry found in Q. 
We would now like to go further and study more generic solutions of the algebraic master equation (13.121) governing 
the general solution of the system. The first important thing to note is that the solution of the cubic k = k{r), as 
given by p.l2l) does not depend on the mass integration constant yx of the solution (13.lip . The mass of the black hole 
solution does not interfere in the equation itself. The algebraic equation (13.121) specifies the asymptotic behavior of 
the solution via the integration constant Cq. As such, solution p.l3p corresponds to de Sitter space with a non trivial 
scalar field </> and can be promoted to a black hole for /x 7 ^ 0. We will refer to it as a de Sitter solution from now on 
(see in that aspect 0). What other solutions fc = fc(r) of de Sitter asymptotics exist other than (13.131) ? It is easy to 
verify that asymptotically (anti) de Sitter solutions will be obtained by fixing the constant Cq to its found de Sitter 
value given above for the de Sitter Schwarzschild geometry (13.131) . The fake hair charge q remains a free parameter. 
We will look into these solutions in order to understand their behaviour beyond the GR stealth limit (I3.13P and in 
particular to see locally around (13.131) how self tuning is influenced in this class of solutions. 

But first of all there is also a second class of solutions we can obtain relatively easily. They are simply given by 
dividing the third order polynomial by the de Sitter root (13.131) . The resulting quotient polynomial will yield a second 
degree equation for fc(r) which can have real roots and which can be easily solved. Indeed we know that the parameter 
choice q^ = and Cq = ( 2 ^ 7 ? — \)\fP/'q yields (13.131) as a solution of (13.121) . 

In what follows we will consider the cases /3 < 0 and /3 > 0 separately. For the second order equation it is /3 > 0 
that presents more interest in terms of regular solutions. For the former de Sitter solutions they will be in the same 
family as (13.131) hence /3 < 0. Now that we have understood the influence of each term of the action in the field 
equations we will set natural values for the couplings in the action = 1, q = 1/2 and consider A > 0 while defining 
R = -ylpj and u = 2|/3|A without loss of generality. Hence the u parameter is dimensionless and denotes the strength 

of the higher order coupling /3 for given cosmological constant. 






































(4.1) 


For f3 > 0 we have q‘^\l3\ = (1 + u), Co = — u) and setting 


fcl/2 


the cubic (13.121) takes the canonical form, 




2 + {l-u)R^ 


(1 - u) 


(l+M)(l + i?2)2 (l+w)(l + i?2)2 

Solution (13.131) reads in the novel rescaled variables, 

1 


= 0 


l + i?2 


So, writing again the cubic as 


we conclude that 


while 


where we have set for simplicity 


~ TT^) + bx + c)=0 


b = 


1 


l + i?2 


1 


l + i?2 


> 0 


= —vb 


V = 


I — u 

1 + M 


It is then easy to show that the remaining two solutions of (14.21) are given by 

1 


a:i,2(.R) = 


2(l + i?2) L 


-l±\/l + 4u(l + i?2) 


(4.2) 


(4.3) 


(4.4) 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


For u > 0 both of the above roots are real which is true for —1 < 2A/3 < 1 and the solution avoids branch singularities 
(where the square root is zero). Computing the integral in p.lll) we get. 


R h{r) = 


k{r) 


P + rjr 


-dr = 4 


(l+i?2) 


-1± V'l + 4i;(l + i?2) 


jdR 


(4.9) 


whereas 


f{R) = h{R) 


1± \/l + 4u(l + i?2) 


2 


(4.10) 


Therefore / has at least the same zeros as the function h which will be event horizons. At the end the solution for h 
reads after direct integration, 

h{R) = i — — H-i— [ arctan(i?) ± arctan ( — ^ ] | ± — sinh”^ ( J ] (4-11) 

^ ’ V R 2v-^R\ ^ ’ \^v/l + 4i’(i?2 + l ))) [yi + AvJ '' ^ 

Therefore, for large R, h = h{R) is bounded and asymptotes 1/v. For A = 0 we have u = 1 and this is the only case 
where the area of the sphere is not reduced. Otherwise we have a solid angle deficit which at the absence of an event 
horizon can lead to a naked singularity. For fi > 0 there is always an event horizon covering the singularity at i? = 0 
and the above solution is a black hole with a well behaved scalar field. We have found therefore a regular scalar tensor 
black hole solution quite similar to the black hole embedded in an Einstein static universe Q. 
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A second case of interest to consider is to take Cq = (2C?7 — X)y/^/ri while allowing for arbitrary charge q. Our aim 
here is to seek solutions which asymptote de Sitter space and generalize the de-Sitter Schwarzschild solution found in 
Q given by (I3.13|) . So in accordance to this aim we take /3 < 0, that is /3 = —1/3|. Now, the constant Cq takes the 
value 


so that 


Co=^^^(l + ^i) 



(4.12) 


(4.13) 


as r —> oo. Given the imaginary value of Cq, setting ix 


VW\ 


the cubic becomes 


2-(l + u)i?2 (l+rt) 


(4.14) 


Expressions are slightly more involved here and we make some redefinitions to render expressions more readable. 
Setting, 


the cubic reads, 


?/ = 1 - 


u + 1 


y + v _^ 

p2y2 ^ 


= 0 


where v is given in (14.711 and we note that, 


where XQ{y) is solution to the cubic (|4.16|) . Customary to third order equations we define 

P={y + ’v)Q 

p y 

The number of real solutions of (14.141) depends upon the sign of the discriminant 


(4.15) 

(4.16) 

(4.17) 

(4.18) 

(4.19) 


More specifically , for A > 0 we have only one real root and a complex pair. Since we are looking for a static patch 
of De Sitter we take 0 < y < 1. Indeed the de Sitter Schwarzschild solution is obtained for = —v and Xds ~ y ■ 
is then clear that A > 0 and the unique real solution is given by 


Xq = 


2^/3py 


-py + yp'^y'^ + ^( 2 / + 


1/3 


2l/3(y + „) 


Spy 


-py 


p2y2 


{y + vf 


1/3 


(4.20) 


The explicit solution can be obtained by a simple but tedious procedure. Here we are mostly interested in deter¬ 
mining the asymptotic behavior of the solution close to the de Sitter like horizon in y = 0. Hence we consider the 
leading order behavior of the real solution(s) in small y for differing values of v respectively negative and positive. 
The stealth de-Sitter solution (13.1311 belongs to the family with u < 0. In this case we obtain, xq ^ and this 
is again a de-Sitter Schwarzschild solution to leading order in y. Therefore, the de Sitter Schwarzschild solution is 
not a particularly special solution it is continuously related to a full branch of de Sitter like black hole solutions with 
similar characteristics. For > 0 we get however xq ^ ^ which yields Rh = tanh”^ R — y which has a singular 
behavior asymptotically. Therefore we see that for u < 0 the solutions are indeed similar to (13.131) . 






























10 


In order to back up this conclusion without going to numerics we can also consider the following approximation. 
Suppose that we have a small deviation from the charge configuration valid around the de-Sitter solution (I3.13L 
namely 

9^I/3| = (m- 1)+ e (4.21) 

with e << |u — 1|. We consider the linearized solution in e to (14.141) . 

X = xo + exi (4.22) 

where xq = "^ 3 ^ and xi is the deviation from the de-Sitter solution (I3.13P when we shift slightly the fake hair 
parameter. Now, upon substituting the above into (14.141) to linear order in e we get, 


x{R) 


1 

(1 - R^) 


(3u- 1) - {l + u)R^ 


Now, if u G (—1,1/3) defining, 


1 -I- u 
1 - 3 m 


7 ^ > 0 


we get 


h{R) = 



Th‘) 


L 

R 


i?2 + 1 

3 (1 -I- u)"f R 


arctan (jR) 


and 


f{R) 


VW\R 



2e \ 2e 

(3m- 1) - (1 -f n)i?2 ) ~ l + u 

(1 -I- m)7 R 


(4.23) 


(4.24) 


(4.25) 


(4.26) 


As long as our approximation e << |m, — 1| is true the linear perturbative solution is valid for large enough R and 
the solution is regular and very similar to a de Sitter Schwarzschild solution. Note that the correction in e does not 
interfere with the effective cosmological constant which means that the self-tuning mechanism is stable under a small 
deviation of the coupling q. Alternatively for u G (—oo, —1) U (1/3, -boo) we can set, 


oj = i"f = 


1 -b M 
3u- 1 


(4.27) 


to get, 


h{r) 



^ R^ 


^ In 

1 -b loR 

R 3 

1 - J 

(1 + M)i?^ 

1-loR 


(4.28) 


where we note now that the solution is again regular and the only difference is the slightly differing behavior of the 
logarithmic term. 


V. CONCLUDING REMARKS 

In this paper we have considered a resolution method originally found in Q giving regular black holes for scalar 
tensor theories. Following the analysis in Q we showed how the complex system of field equations boiled down to 
solving a single third order algebraic equation (I3.12|) which basically fixed the asymptotic properties of the spacetime 
solution. We then presented some new solutions analyzing in some detail the algebraic equation at hand (13.121) which 
depends on two integration constants Cq and q (but not the mass of the black hole). The two integration constants 
are the fake hair charge q and an integration constant Cq determining the asymptotic properties of the solution at 
hand. Apart from novel regular black holes we also found that the particular de Sitter self tuning solution p.l3l) 
found in [8| which is a self tuning black hole solution, is not some isolated point in the space of algebraic solutions, but 
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is a generic and stable configuration with respect to small deviations in the fixed parameters of the solution. Indeed 
when we consider a small deviation in the self-tuning integration constant q the effective cosmological constant of the 
neighboring black hole solution remains the same not spoiling the self tuning mechanism. 

The regularity and integrability of the method presented relies on three properties for the Galileon scalar field: first 
of all it has to have translational invariance and this symmetry yields a conserved charge and current with a special 
form for the scalar field equations. Secondly, the scalar tensor action considered has to include higher order Galileon 
terms so that the scalar held equation admits novel, geometric branches of solutions apart from the ones that lead to 
trivial or singular solutions. This enables us to satisfy the no hair arguments present for Galileons while hnding 
non trivial scalar tensor black holes (see also 0). The regularity requirements required in are still not sufficient 
to obtain regular solutions. We need a third ingredient property : the key to obtaining a regular scalar held at the 
horizon is to additionally allow a linear time dependence for the scalar held. Under these three important properties 
it is then relatively easy to hnd novel scalar tensor black hole solutions. It would be interesting to see how far can one 
generalize these principles to go to more general Horndeski theories (see 4^ for progress in this direction) or again 
hnd hairy black holes with a genuine scalar charge. In the latter case hairy black holes with genuine scalar charge 
have been found in a bi scalar tensor theory that are seen to elegantly generalize the BBMB black hole. Another 
interesting open question lies in the stability of such scalar-tensor black holes as those discussed here. 
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